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ABSTRACT 

In recent years, the evidence for the existence of an ultra-compact concentration of dark mass associated 
with the radio source Sgr A* in the Galactic Center has become very strong. However, an unambiguous 
proof that this object is indeed a black hole is still lacking. A defining characteristic of a black hole 
is the event horizon. To a distant observer, the event horizon casts a relatively large "shadow" with 
an apparent diameter of ~ 10 gravitational radii due to bending of light by the black hole, nearly 
independent of the black hole spin or orientation. The predicted size (~ 30/iarcseconds) of this shadow 
for Sgr A* approaches the resolution of current radio-interferometers. If the black hole is maximally 
spinning and viewed edge-on, then the shadow will be offset by ~ 8 /iarcseconds from the center of 
mass, and will be slightly flattened on one side. Taking into account scatter-broadening of the image in 
the interstellar medium and the finite achievable telescope resolution, we show that the shadow of Sgr 
A* may be observable with very long-baseline interferometry at sub-millimeter wavelengths, assuming 
that the accretion flow is optically thin in this region of the spectrum. Hence, there exists a realistic 
expectation of imaging the event horizon of a black hole within the next few years. 

Subject headings: black hole physics — relativity — Galaxy: center — galaxies: active — 
submillimeter — techniques: interferometric 



1. INTRODUCTION 



High resolution spectroscopy (especially with the Hub- 
ble Space Telescope) of galactic nuclei has produced an 
abundance of evidence for compact dark mass concentra- 
tions of up to 10 8 Mq/pc 3 , whose nature is strongly sus- 



pected to be indicative of supermassive black holes (Ko 



rmendy fc Richstonc 1005 ): — Even better evidence cxiat3 



A* itself, since the latter, unlike the surrounding stars, has 
a tightl y restricted prop e r motion indicating that it is very 
heavy ( |Reid et al. 1000| , |Genzel et al. 1997| ). 

The key spectral features of Sgr A* are a slightly in- 
verted cm-wavelength spectrum, an apparent excess (or 
bump) at sub-millimeter (sub-m m) wavelengths, and a 
steep cut-off towa rds the infrared ( Falcke et al. 1008] , her- 



i tor the galaxy NGC 4258 and the Milky Way, for which 



abyn et al. 1007). The radio emission is cir cularly polar 



spectroscopic and proper motion studies have provided an 
unprecedented throe dimonoional view of the kinomatico of 



ized but undetected in linea r polarization (Bower et al 
1000a , Bower et al. 1000b| ). Proposed models f or the 



ga3 and stars around a central point ma33, pointing to dark 
, I macs conc entrations of > 10 12 J \ / 1 i ) /pc 3 with very high -sig 



radio emission range from quasi-spherical inflows ( J Melia 
1002| , |Melia 1994 iNarayan, Yi, fc Mahadevan lflQ sj) to a 



nificancc ( i Miyoahi ct al. 1995 | , Eckart fc Ccnzel VMdt -ht 



jet-like outflow (Falcke, Mannheim, fc Biermann 1993| , |Fal- 
cke fc Biermann 1999] ) 



ot al. 100 



Complementary observations of galactic nuclei with very 
long baseline interfero metry (VLBI) reveal the presence of 
compact radio cores (Zensus 1007) which appear to be 
coincident with the central black hole candidates. An 
intriguing case is that of the Galactic Center where the 
bright, co mpact radio source S g r A* lies at the dy nami- 
cal origin ( |Menten et al. 1997| , |Ghez et al. 1998^ . The 
nature of Sgr A* is still unclear, since its structure is com- 
pletely washe d out by strong i nterstellar scattering at cm- 
wavelengths ( [Lo ct al. 1008 ). It is only at millimeter- 
wavel cngths that we may begin to sec some internal struc- 
ture (Bower fc Backer 1Q08| , Krichbaum ct al. 1008 , Lo ct 



al. 1008). Though the dark mass concentration could in 



principle be distributed in the form of exotic objects on a 
scale slightly larger than the size of Sgr A* — but wi th dif- 
ficulties accoun ting for its radiation characteristics ( Melia| 
fc Coker 1000) — it is expected to be associated with Sgr 



The sub-mm bump is particularly interesting since this 
should be the signature of a very compact synchrotron 
emitting region with a size of a few Schwarzschild radii 
flFalckc 1996[ [Falcke ct al. 1998| ). The presence of com- 
pact radio emission in Sgr A* at a wavelength as short as 
1.4 mm has been confirme d recently by a first VL BI de- 
tection at this wavelength ( Krichbaum et al. 1995 ). This 
detection is exciting for several reasons. First, it lies in 
a region of the spectrum where the intrinsic source size 
should become ap parent over scatter-broadening by the in- 
tervening screen (Melia, Jokipii, & Narayanan 1902). Sec- 
ond, this component is sufficiently bright to be detected 
with VLBI techniques at even shorter wavelengths, and 
third, Sgr A* is sufficiently close that the size scale where 
general relativistic effects are significant could be resolved 
with VLBI at sub-mm wavelengths. In addition, at sub- 
mm wavelengths, the various models predict that the syn- 
chrotron emission is not self-absorbed, allowing a view into 
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the region near the horizon. The horizon has a size of 
v/l — al)R g , where R g = GAI/c 2 , M is the mass of 



1 



the black hole, G is Newton's constant, c the speed of light, 
a* = Jcj (GM 2 ) is the dimensionless spin of the black hole 
in the range to 1, and J is the angular momentum of the 
black hole. 

Bardeen (1973) described the idealized appearance of a 
black hole in front of a planar emitting source, showing 
that it literally would appear as a 'black hole'. At that 
time such a calculation was of mere theoretical interest 
and limited to just calculating the envelope of the appar- 
ent black hole. To test whether the re is a realistic ch ance 
of seeing this 'black hole' in Sgr A* ( Falcke et al. 1998 ), we 



here report the first calculations obtained with our general 
relativistic (GR) ray-tracing code that allows us to simu- 
late observed images of Sgr A* for various combinations of 
black hole spin, inclination angle, and morphology of the 
emission region directly surrounding the black hole and not 
just for a background source. A mor e detailed description 
of our calculations is in preparation ( Agol, Falcke, fc Melia 
19990 

2. THE APPEARANCE OF A BLACK HOLE 

We determine the appearance of the emitting region 
around a black hole under the condition that it is optically 
thin. F or Sgr A* this mig ht be the case for the sub-mm 
bump ( Falcke et al. 199q ) indicated by the turnover in 
the spectrum, and can always be achieved by going to a 
suitably high frequency Here we simply assume that the 
overa ll specific intensity, observed at infinity is an in- 



tegrat ion of the cmissivity, j v , times the differential path 
lengt h along geodesies ( |Jaroszynski fc Kurpiewski 1997|) . 
In line with the qualitative discussion of this paper, we 
assume that j u is independent of frequency, and that it is 
either spatially uniform, or scales as r~ 2 . These two cases 
cover a large range of conditions expected under several 
reasonable scenarios, be it a quasi-spherical infall, a rotat- 
ing thick disk, or the base of an outflow. 

The calculation of the photon trajectories and the in- 
tensity integrated al ong the linc-o f -sight is based on the 
standard formalism ( Thome 1981 j_ Vicrgutz 1993 j_ Rauch| 
fe Blandford 1994 paroszynski fc Kurpiewski 1997 ). Our 
calculations take into account all the well-known relativis- 
tic effects, e.g.. frame dragging, gravitational redshift, light 
bending, and Doppler boosting. The code is valid for all 
possible spins of the black hole and for any arbitrary ve- 
locity field of the emission region. 

For a planar emitting source behind a black hole, a 
closed curve on the sky plane divides a region where 
geodesies intersect the hor izon from a region whose 
geodesies miss the horizon ( Bardeen 1973 ). This curve, 
which we refer to as the "apparent boundary" of the black 



hole, is a circle of radius v27R g in the Schwarzschild case 
(a* = 0), but has a more flattened shape of similar size for 
a Kerr black hole, slightly dependent on inclination. The 
size of the apparent boundary is much larger than the 
event horizon due to strong bending of light by the black 
hole. When the emission occurs in an optically thin re- 
gion surrounding the black hole, the case of interest here, 
the apparent boundary has the same exact shape since 
the properties of the geodesies are independent of where 
the sources are located. However, photons on geodesies 



located within the apparent boundary that can still es- 
cape to the observer experience strong gravitational red- 
shift and a shorter total path length, leading to a smaller 
integrated emissivity, while photons just outside the ap- 
parent boundary can orbit the black hole near the circular 
phot on radius several times, adding to the observed inten- 
sity ( Jaroszynski fe Kurpiewski 1997 ). This produces a 
marked deficit of the observed intensity inside the appar- 
ent boundary, which we refer to as the "shadow" of the 
black hole. 

We here consider a compact, optically-thin emitting re- 
gion surrounding a black hole with spin parameter a* = 
(i.e., a Schwarzschild black hole) and a maximally spin- 
ning Kerr hole with a* = 0.998. In the set of simulations 
shown here, we take the viewing angle i to be 45° with 
respect to the spin axis (when it is present), and we con- 
sider two distributions of gas velocity v. The first has 
the plasma in free-fall, i.e., v r — — ^/2r(a 2 + r 2 )A/A and 
2ar/A, where v r is the Boyer-Lindquist radial ve- 



locity, is the orbital frequency, A = r 2 — 2r 



and 



A = (r 2 + a 2 ) 2 - a 2 A sin 9. (We have set G = M = c = 1 
in this paragraph.) The second has the plasma orbiting 
in rigidly rotating shells with the equatorial Keplerian fre- 
quency f2 = l/(r 3 / 2 + a) for r > r ms with v r = 0, and 
i nfalling with const ant angular momentum inside r < r ms 
( Cunningham 1975|) , with v e = for all r. 

In order to display concrete examples of how realistic 
our proposed measurements of these effects with VLBI will 
be, we have simulated the expected images for the mas- 
sive black hole candida te Sgr A* at the Ga lactic Center. 



For it s measured mass ( Eckart fc Genzel 1996| , |Ghcz et al 
1998D M = 2.6 x 10 6 M , the scale size for this object is 



the gravitational radius R g = 3.9 x 10 11 cm, which is half 
of the Schwarzschild radius R s = 2GM/c 2 . 

To simulate an observed image we have to take two ad- 
ditional effects into account: interstellar scattering and 
the finite telescope resolution achievable from the ground. 
Scatter-broadening at the Galactic Center is incorporated 
by smoothing the image with an elliptical Gaussian with a 
FWHM of 24.2 /iarcsecondx (A/1.3 mm) 2 along the major 
axis and 12.8 /iarcse condx (A/1.3 mm) 2 along the minor 
axis ( Lo et al. 1998| ). The position angle of this ellipse 
is arbitrary since we do not know yet the spin axis of the 
black hole on the sky and we have assumed PA=90° for 
the major axis. The telescope resolution — in an idealized 
form — is then added by convolving the smoothed image 
with a spherical Gaussian point-spread function of FWHM 
33.5 ^arcsecondx(A/1.3mm) -1 (^/8000km) -1 — the possi- 
ble resolu tion of a global interferometer with 8000 km 
baselines ( Krichbaum 1996 ). In reality the exact point- 
spread-function will of course depend on the number and 
placement of the participating telescopes. 

In Figure 1, we show the resulting image of Sgr A* 
for a maximally rotating black hole viewed at an angle 
of i = 45°, for a compact region in free fall, with an emis- 
sivity j u — is°r~ 2 . We first show the original, unsmoothed 
image of the emission region as calculated with the GR 
code in panel (a), and then present the simulated 'ob- 
served' images at 0.6 and 1.3 mm wavelengths in panels 
(b) and (c), respectively. The two distinct features that 
are evident in Figure la are (1) the clear depression in I v — 
the shadow — produced near the black hole, which in this 
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Fig. 1. — An image of an optically thin emission region surrounding a black hole with the characteristics of Sgr A* at the Galactic Center. 
The black hole is here either maximally rotating (a* = 0.998, Figs, la-c) or non-rotating (a* = 0, Figs. ld-f). The emitting gas is assumed to 
be in free fall with an emissivity oc r -2 (top) or on Keplerian shells (bottom) with a uniform emissivity (viewing angle i = 45°). Figs, la&d 
show the GR ray-tracing calculations, Figs, lb&e are the images seen by an idealized VLBI array at 0.6 mm wavelength taking interstellar 
scattering into account, and Figs, lc&f are those for a wavelength of 1.3 mm. The intensity variations along the x-axis (solid green curve) 
and the j/-axis (dashed purple curve) are overlayed. The vertical axes show the intensity of the curves in arbitrary units and the horizontal 



axes shows the distance from the black hole in units of R s which for Sgr A* is 3.9 X 10 



3 ^arcseconds. 



particular example represents a modulation of up to 90% 
in intensity from peak to trough, and (2) the size of the 
shadow, which here is 9.2i? g in diameter. This represents a 
projected size of 27 /xarcseconds, which is al ready within a 
factor of two of the current VLBI resolution ( Krichbaum et 
al. 19D3). The shadow is a generic feature of various other 
models we have looked at, including those with outflows, 
cylindrical emissivity, and various inclinations or spins. 

To illustrate the expected image for another extreme 
case, we show in Figure Id the analogue to Figure la for 
the case with a* = (i.e., no rotation), an emitting plasma 
orbiting in Keplerian shells (as described above), and a 
uniform j v for r < 25 R g . Even though these conditions 
are distinctly different compared to those of Figure la, the 
black hole shadow is still clearly evident, here representing 
a modulation in I u in the range of 50-75% from peak to 
trough (Fig. Id), and with a diameter of roughly 10.4 R g . 
In this case, the emission is asymmetric due to the strong 
Doppler shifts associated with the emission by a rapidly 
moving plasma along the line-of-sight (with velocity v$). 

The important conclusion is that the diameter of the 
shadow — in marked contrast to the event horizon — is fairly 
independent of the black hole spin and is always of or- 
der 10i? g . Indeed, this is consistent with the observed 
0.8 mm size limit > AR a of Sgr A* from a lack of scin- 
tillation (Gwinn et al. 1991). The presence of a rotating 
hole viewed edge-on will lead to a shifting of the apparent 
boundary (by as much as 2.5 R g , or 8 /xarcseconds) with 
respect to the center of mass, or the centroid of the outer 



emission region. 

Interestingly, the scattering size of Sgr A* and the res- 
olution of global VLBI arrays become comparable to the 
size of the shadow at a wavelength of about 1.3 mm. As 
one can see from Figures lc&f the shadow is still almost 
completely washed out for VLBI observations at 1.3 mm, 
while it is very apparent at a factor two shorter wave- 
length (Figures lb&e). In fact, already at 0.8 mm (not 
shown here) the shadow can be easily seen. Under certain 
conditions, i.e., a very homogeneous emission region, the 
shadow would be visible even at 1.3 mm (Fig. If). 

3. HOW REALISTIC IS SUCH AN EXPERIMENT? 

The arguments for the feasibility of such an experiment 
are rather compelling. First of all, the mass of Sgr A* is 
very well known within 20%, the main uncertainty be- 
ing the exact distance to the Galactic Center. Since, 
as we have shown, the unknown spin of the suspected 
black hole contributes only another 10% uncertainty, we 
can conservatively predict the angular diameter of the 
shadow in Sgr A* from the GR calculations alone to be 
~ 30 ± 7^arcseconds independent of wavelength. As seen 
in Fig. 1, the finite telescope resolution and the scatter 
broadening will make the detectability of the shadow a 
function of wavelength and emissivity; however, the size 
of the shadow will remain of similar order and under no 
circumstances can become smaller. 

The technical methods to achieve such a resolution at 
wavelengths shortwards of 1.3 mm are currently being de- 
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veloped and a first detection of Sgr A* at 1.4 mm with 
VLBI has already been reported. The challenge will be to 
push this technology even further towards 0.8 or even 0.6 
mm VLBI. Over the next decade many more telescopes 
are expected to operate at these wavelengths. Depend- 
ing on how short a wavelength is required, the projected 
time scale for developing the necessary VLBI techniques 
may be about ten years. A fundamental problem prevent- 
ing such an experiment is not now apparent, but in light of 
our results, planning of the new sub-mm-telescopes should 
include sufficient provisions for VLBI experiments. 

A potential problem with our model may occur if j u 
has an inner cutoff which is larger than that of the hori- 
zon, making the shadow larger than predicted due to a 
decrease in emissivity rather than to GR effects. How- 
ever, first of all, the truncation of accretion disk emis- 
sion a t the marginal stab le orbit r ms is somewhat arbi- 



trary ( Cunningham 1975| ) and, secondly, if it exists such 
a cutoff would likely be frequency dependent, while there 
will be a frequency-independent minimum radius due to 
the general relativistic effects we have described. Another 
problem could be the unknown morphology of the emis- 
sion region. Anisotropy, strong velocity fields, and den- 
sity inhomogeneities would make an identification of the 
shadow in an observed image more difficult. However, in- 
homogeneities are unlikely to be a major issue, since the 
time scale for rotation around the black hole in the Galac- 
tic Center is only a few hundred seconds and hence much 
less than the typical duration of a VLBI observation. The 
strong shear near the black hole would tend to smooth 
out any inhomogeneities very quickly. Indeed, sub-mm 



variability studies on such short time scales ( Gwinn et al_ 
199l[j]have yielded negative results. The same argument 
applies to emission models which are offset from the black 
hole, e.g., are one-sided. Since the shadow of the black hole 
has a very well defined shape it would under any condi- 
tions appear as a distinct feature, given that the dynamic 
range of the map is large enough (i.e., J>100:1, considering 



a range of emission models, Agol, Falcke, & Melia 1999). 

Finally, synchrotron self-absorption could pose a prob- 
lem. So far the available sub-mm spectra show a flattening 
of the spectrum around 1.3-0.6 mm indicating a turnover 
towards an optically thin spectrum. Given the current ob- 



servational uncertainties one could in principle construct 
simple models where the flow does not become optically 
thin until 0.2 mm. Improved simultaneous measurements 
at sub-mm wavelengths are therefore highly desirable to 
exactly measure the spectral turnover since the experi- 
ment wc propose here will only work for an optically thin 
flow. At hundreds of microns the atmosphere becomes op- 
tically thick, making much more expensive space-based ob- 
servations necessary. At X-ray wavelengths, the accretion 
flow will be optically thin to electron scattering, so there 
may be a better chance of detecting the shadow with fu- 
ture space-based X-ray interferometry as proposed in the 
MAXIM experiment. 

4. SUMMARY 

The importance of the proposed imaging of Sgr A* at 
sub-mm wavelengths with VLBI cannot be overempha- 
sized. The bump in the spectrum of Sgr A* strongly sug- 
gests the presence of a compact component whose proxim- 
ity to the event horizon is predicted to result in a shadow 
of measurable dimensions in the intensity map. To our 
knowledge, such a feature is unique and Sgr A* seems 
to have all the right parameters to make it observable. 
The observation of this shadow would confirm the widely 
held belief that most of the dark mass concentration in 
the nuclei of galaxies such as ours is contained within 
a black hole, and it would be the first direct evidence 
of the existence of an event horizon. A non-detection 
with sufficiently developed techniques, on the other hand, 
might pose a major problem for the standard black hole 
paradigm. Because of this fundamental importance, the 
experiment we propose here should be a major motiva- 
tion for intensifying the current development of sub-mm 
astronomy in general and mm- and sub-mm VLBI in par- 
ticular. 
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